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Our peculiar motion with respect to the CMB rest frame represents a preferred direction in the 
observed CMB sky since it induces an apparent deflection of the observed CMB photons (aberration) 
and a shift in their frequency (Doppler). Both effects distort the multipoles a^m's at all ^'s. Such 
effects are real as it has been recently measured for the first time by [1] according to what was 
proposed in [2-4] . However, the common lore when estimating a power spectrum from CMB is to 
consider that Doppler affects only the £ = 1 multipole, neglecting any other corrections. In this letter 
we use simulations of the CMB sky in a boosted frame with a peculiar velocity /3 = v/c = 1.23 x 10"'' 
in order to assess the impact of such effect on power spectrum estimations in different regions of 
the sky. We show that the boost induces a north-south asymmetry in the power spectrum which 
is highly significant and non-negligible, of about (0.75 ± 0.15)% for half-sky cuts when going up to 
£ ~ 2500. We suggest that these effects are relevant and may account for some of the north-south 
asymmetries seen in the Planck data, being especially important at small scales. 

Keywords: CMB theory, CMB aberration, CMB anomalies 



o 
u 

Oh 
I 

o 



> 

o 
in 

en 

^' 

o 
en 

> 

X 



Introduction. The Cosmic Microwave Background is 
used as a fundamental tool to test Cosmological Mod- 
els and to quantitatively estimate the parameters of such 
models. This is usually done by extracting the Temper- 
ature and the Polarization power spectra from the maps 
and by fitting them with a ACDM background model 
with an almost scale-invariant gaussian spectrum of den- 
sity perturbations. However, if we observe the CMB with 
a velocity /? = v/c relative to such background, the image 
undergoes distortions due to the Doppler and aberration 
effects and the correct procedure should be to first trans- 
form the image in the CMB rest frame and then analyze 
the data. Surprisingly such procedure is not performed 
by any of the current experimental analysis. 

Such boost in a direction defined by the unit vector z 
distorts a primordial temperature map r(n) by a Doppler 
factor and by changing the arrival direction: 



T'(n')-7(l + /3n-z)T(n), 



(1) 



where (fi' — n)-z 



Psin^ 6/{l+Pcos9), cos6 = n-z, and 
7 = i/-\/i — /32. Such distortion introduces correlations 
between different a^^'s (the coefficients of the spherical 
harmonics decomposition), both diagonally (i.e., between 
same ^'s) and ojf-diagonally (between different ^'s) [5]: 
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The only effect which is usually taken into account is 
the large £ — 1 dipole, due to the matrix element Kqiq, 
which is used to infer our velocity /3 [6]. Clearly how- 
ever there is much more information in the other ma- 
trix elements. The off-diagonal correlation was shown 
in [2-4] to be measurable by Planck as an alternative 
method to measure f3 and in fact the Planck collabora- 



tion itself [1] recently published for the first time a de- 
tection of /3 through CMB aberration and Doppler using 
1 < £ < 2000. This provides an alternative measurement 
of /3 with respect to the usual £ = 1 Doppler effect. For 
the diagonal correlation it has been shown in a previous 
paper ]7] (see also ]8] for earlier work) with numerical 
simulations of the CMB sky that while in the full-sky 
case it is safe to disregard such boost effects, there is 
instead a non-negligible distortion in power spectrum re- 
construction due to the presence of a mask which may 
induce a significant bias in the parameter estimation for 
Planck, up to 1 — 2 standard deviations. 

In this letter we assess whether ignoring such frame 
effect may also induce a significant spurious directional 
dependence in the CMB power spectrum. In particular 
we analyze the pseudo-C^ reconstructed in different re- 
gions of the sky by performing simulations in the case 
in which a boost effect is present or not. We proceed 
here as in ]7]: we simulate maps of the CMB sky and 
we directly apply on the maps the boost transformation 
before extracting the aim's, therefore bypassing the need 
of computing the mixing coefficients. Then we extract 
the agm's and the C'e's and show that an asymmetry is 
visible. We also roughly quantify the amount of asym- 
metry through an asymmetry coefficient which is built to 
give an estimate of the change of the overall amplitude 
of the power spectrum. We apply all this procedure to 
Temperature maps for an ideal experiment (neglecting in- 
strumental noise) for 3 different symmetric sky-cuts with 
measurable fraction of the sky /sky = {0.88, 0.29, 0.04}) 
to a maximum multipole £max = 2700. Our results can 
thus easily represent a Planck-like case (almost all-sky 
and ^max — 2500), a WMAP-like case (almost all-sky 
and £max — 900) and in part some other current small- 
scale surveys, like SPT [9] and ACT [10] (smaller /^ky 



and Cax=i 3000). 

Power spectra asymmetry. We apply our procedure to 
extract pseudo-C^'s from different regions of the sky in 
simulations performed with the HEALPix package^ in a 
modified version [7] which allows boosts, with A'sidc — 
4096 and ^max = 2700. Such a numerical boost proce- 
dure has been tested with Bessel fitting functions which 
reproduce the K^iirn with high precision [4, 7\? We ap- 
ply a boost along the north pole direction of /3 = 0.00123 
and compare two regions in two opposite directions along 
such boost (dubbed North and South). We performed 
20 different simulations as random realizations, labelled 
by the random seed used as an input in HEALPix, of a 
fiducial cosmological model similar to the WMAP 9-year 
best fit. We find that there is a systematic asymme- 
try between power in the two opposite directions in the 
case with boost as opposed to the case without boost. 
We show the difference in power spectra in Fig. 1. This 
should be compared with the experimental results pre- 
sented in [12] (section 5.5.1 and Fig. 28), from which 
we can see that the real data has systematically more 
power in one hemisphere than the other at a level of a 
few percent, although it is difficult to make a more pre- 
cise estimate on the significance and scale of the effect 
based on the results presented at that paper. 

Given a set of C^'s it is possible to have a rough esti- 
mate on the size of the effect on a cosmological param- 
eter, by considering [7] an idealized case in which the 
CMB depends multiplicatively on a single amplitude pa- 
rameter, which we call A, so that the y^ is given by: 
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where Cf^ is the theoretical spectrum when A—1 
are the observed values in one region of the sky and where 
of = C| 2/(2£ -t- 1) is the cosmic variance, ignoring any 
noise. The best fit value Abf for A is obtained when 

.2\ 



d{x )/dA = which gives: 



A 



bf 



E 






E 



{C 



th\2 



(4) 



The difference A^ 



\i- 



between the best-fit values in 
two regions N and S of the sky with observed spectra 
C™''' and C^'^^' is therefore given by: 



5_A 
~A 



^bf + ^bf 



,?f'^'-> 




^ http : //healpix . sourcef orge . net/ 

^ We have checked up to £max = 3000 that boosting a A^sidc = 4096 
map with /3 = 0.00123 and boosting again with — /3 gives back 
the original map with very high precision up to the 7*'^ digit in 
the a^^'s, differently from what obtained in [11]. 
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Figure 1. Relative difference between Ce's in two opposite 
discs of the sky centered on the dipole direction as a function 
of the multipole £. The brown (green) line shows the mean 
spectrum over 20 boosted (unboosted) simulations, binned in 
100-^ bins. The green area shows the la band around the 
unboosted mean. [Top]: two /salves of the sky (/sky ~ 0.44, 
after removing a band around the galaxy). [Middle]: two 
antipodal discs of 90° diameter (/sky = 0.146, to be compared 
with [12]). [Bottom]: 4 random boosted realizations. 



where 6Ce = C^^^' — C™'^' and we have approximated 
C|'' by the average (Cf ^'^ -I- C°''^'^)/2. We therefore 
estimate^ SA/A for our simulations summing up to a cer- 
tain (-max- We will compute the estimator directly us- 
ing the pseudo-Cf's given as an output by HEALPix in 
eq.(5). For a sufficiently large patch of the sky the off- 
set between pseudo-C^ and Ci is simply an overall factor 
given by the sky fraction /sky so our eq.(5) is still a good 
estimator of the difference in amplitude of the two hemi- 
spheres. Note however that we will also use the same 



^ Such an estimator has been used also by [13] to quantify hemi- 
spherical asymmetries. 
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Figure 2. Average and la bands for 20 realizations of the 
North-South asymmetry parameter SA/A of eq.(5). Brown 
curves are made on boosted maps, green on non-boosted ones. 
The difference is computed using pseudo-C^'s in two halves 
of the sky (/sky = 0.44, after removing a band around the 
galaxy). /Tbp/ boost towards the north pole. [Bottom]: boost 
along the actual direction given by the measured £ = 1 CMB 
dipole, [{l,b) — (264,48)] in galactic coordinates. 



equation for smaller sky patches with the caveat that in 
this case SA/A is not directly related to the change in 
amplitude (since the pseudo-C^ and the Ci are related 
by a window function [14] and also eq.(3) should have a 
full covariance matrix) but still represents a simple mean- 
ingful way to quantify an asymmetry. The typical values 
we obtain for SA/A are small and centered around zero 
for the case without aberration while they represent an 
important effect of order 1% when aberration is present. 
We show in Figs. 2 and 3 the result for mean and stan- 
dard deviation of the asymmetry parameter as a function 
of ^max for our 20 simulations. Such an asymmetry is def- 
initely detectable: at imax ~ 2000, it amounts to 3ct for 
/sky = 2 X 0.44 and to over 3.5a for /g^y = 2 x 0.146, and 
the significance increase with ^max (e.g. it is above 5a at 
^max ~ 2500) . Actually such an asymmetry in the power 
spectrum constitutes by itself a self-consistency check on 
the measurement of /3 given in ]1], as already proposed 
initially by [15]. Given such a large result it is thus natu- 
ral to wonder whether this could account for the asymme- 
tries measured by Planck at 2 < £ < 1500, or at least for a 



Figure 3. Same as Fig. 2 for smaller patches of the sky. Note 
that although the discrepancy is larger for smaller /sky, this 
is over-compensated by the increase of variance among dif- 
ferent realizations. Nevertheless the observed CMB is one 
particular realization and so the difference on a parameter es- 
timation may well be large (of order a few percent). [Top]: 
two antipodal discs of 90° diameter (/sky = 0.146, similarly to 
the very recent Planck collaboration asymmetry analysis (Fig. 
28 in [12]). ]Bottom]: two antipodal discs of 32.5° diameter 
(/sky = 0.02, similar to the partial SPT survey). 



substantial fraction of it. Looking more closely at the re- 
sults of ]12], sect. 5.5.1 we can see that Planck detects an 
overall preferred direction with a significance quantified 
by the fraction of simulations that have smaller cluster- 
ing of the dipole directions than the data: such number 
is of about 11/500, 4/500 or less than 1/500, depending 
on the foreground cleaning methods, which correponds 
to a 97% - 99% C.L. When looking at our Figs. 2 and 3 
at £ = 1500 we can see that a 1.5 — 2cr effect could be eas- 
ily achieved and it may therefore happen that, once the 
boost is subtracted, the Planck anomalies could become 
much less significant or in any case considerably change. 
There is some small effect already at imax — 600 — 900, 
which should be compared with the asymmetries seen 
in the WMAP experiment at 2 < £ < 600 by ] 16- 18], 
see also ]19]. In this case the claimed significance was 
of about 99%C.L. [16-18[ and again it is important re- 
assess such significance once the boost is subtracted. It 
is also interesting to note that the preferred direction 



of WMAP in [16-18] does not quite point in the same 
direction as the dipole but it is not very far way from 
it. More interestingly, for Planck the overall direction 
of the asymmetry is still close to WMAP but it is a bit 
shifted towards the dipole direction. In fact, when look- 
ing at Fig. 27 of [12], which shows the preferred direc- 
tion when analyzing the data in 100 — (. bins, we can 
see that most bins, especially at low i, point towards the 
same direction of WMAP but there are several bins which 
point instead to the dipole direction, especially at large 
L This seems consistent with our findings that at least 
a sizable fraction of the asymmetry may be due to our 
motion in the direction of the dipole, especially at large 
L Note instead that asymmetries and anomalies at very 
large scales ^ ^ 60 [16, 17] cannot be accounted for by 
a boost alone. It is reasonable therefore to consider the 
possibility that there may be some intrinsic large scale 
asymmetry pointing to a different direction which adds 
up with our boost effect. This hypothesis deserves a more 
thorough analysis of the real data, going also to £ > 1500, 
and crucially stresses the need of analyzing the WMAP 
and especially Planck data by first removing the aberra- 
tion and Doppler effects (from the whole map, not only 
the dipole) and then looking for the true significance of 
eventual residual anomalies. 

Fig. 2 (bottom panel) shows the result in the case in 
which the hemispheres are not aligned with the dipole. 
This was achieved by boosting the map along the North- 
South axis and then rotating it with HEALPix along the 
actual dipole direction [{I, h) — (264, 48)] in galactic coor- 
dinates] as measured by the 1=1 multipole of the CMB 
by WMAP and then performing the cut around the North 
and South directions. The results we find for the asym- 
metry parameter 5 A/ A is of lower size and significance, 
as is to be expected since the North-South direction now 
is not anymore the one which maximizes the asymme- 
try. Finally Fig. 3 focuses on the case of smaller discs, 
the first of which allows a direct comparison with ]12] (it 
corresponds to two antipodal 90°-diameter discs). Note 
that in these latter cases the size of the effect is larger, 
because we are selecting the regions of the sky most af- 
fected by a boost (small discs around the poles), but 
at the same time the variance is bigger because we are 
sampling a small fraction of the sky. The bottom panel 
depicts small area analysis that would be relevant for ex- 
periments such as SPT ]9] and ACT ]10] which look at 
small regions (although their actual observation area are 
not circular nor antipodal) of the sky and are therefore 
subject to a potentially large bias due to the boost. We 
leave the analysis of such bias for future work. 

Conclusions. In this letter we have considered the ef- 
fect of our peculiar velocity on the CMB power spec- 
tra in regions in opposite directions of the sky. As in 
]7], we have applied a boost transformation (1) to this 
frame directly in pixel space on simulated maps, rather 
than on the a^m's, which yields a direct test (although 



a similar result could be obtained with the fitting func- 
tions in [4]). We find that aberration and Doppler ef- 
fects induce a directional dependence of the amplitude of 
the power spectrum in opposite directions which we esti- 
mate to be equal to (0.75 ±0.15)% for half sky cuts when 
summing up to ^ < 2500 and is therefore highly signif- 
icant. We claim this could be important to understand 
the Hemispherical asymmetry found by Planck ]12]; at 
least it must account for an important fraction of it. It 
may have also some impact on previous detections in the 
WMAP data at 2 < £ < 600 [16-18], but by itself it 
cannot explain an anomaly localized only at low-£ (say, 
£ < 60 as in ]20, 21]). Note that the direction of max- 
imal asymmetry detected in [12] and [16-18] does not 
coincide with the dipole direction, but some of the bins 
used for the analysis in [12] do indeed point towards the 
dipole direction: this suggests that the asymmetries re- 
ported by Planck are likely to be a mixture of a boost 
effect, showing up especially at high £, with some residual 
asymmetry mostly relevant at low i. It would be useful 
in this respect to analyze the real Planck data at higher 
£ to check whether the asymmetry converges towards the 
dipole direction. Note that we are relying on the stan- 
dard assumption that our velocity is given by the CMB 
dipole; if that is not the case there might be a discrep- 
ancy between the dipole and aberrated directions. In any 
case our findings clearly indicate that before looking for 
preferred directions in the real data at high multipoles 
one should properly deboost the map going to the CMB 
rest frame to avoid spurious detection of anomalies. 

Some of the results in this paper have been derived us- 
ing the HEALPix package ]22]. We thank Jordi Miralda 
and Jorge Noreha for useful discussions and comments. 
MQ is grateful to Brazilian research agencies CNPq and 
PAPER J for support. AN is supported by the grants EC 
FPA2010-20807-C02-02, AGAUR 2009-SGR-168. 
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